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A recently introduced class of quantum spherical spin models is considered in detail. Since the 
spherical constraint already contains a kinetic part, the Hamiltonian need not have kinetic term. As 
a consequence, situations with or without momenta in the Hamiltonian can be described, which may 
lead to different symmetry classes. Two models that show this difference are analyzed. Both models 
are exactly solvable and their phase diagram is analyzed. A transversal external field leads to a 
phase transition line that ends in a quantum critical point. The two considered symmetries of the 
Hamiltonian considered give different critical phenomena in the quantum critical region. The model 
with momenta is argued to be analog to the large- A/" limit of an SU(7V) Heisenberg ferromagnet, and 
, the model without momenta shares the critical phenomena of an SU(A/') Heisenberg antiferromagnet. 

\0 ■ I- INTRODUCTION 

The classical spherical model was conceived by Kac. After being introduced in 1947 to Onsager's rather intricate 
jj» ', solution of the 2d-Ising model, he desired to formulate a simpler spin model. As a first step he took the spins 
' to be continuous Gaussian variables, nowadays called the Gaussian model. This had unphysical behavior at low 
[ temperatures which led Kac to consider the "spherical model" . The spherical model has continuous spins that are 
restricted by the "spherical" constraint X^i^i — which represents the hypersphere intersecting all vertices of 
' the hypercube sustained by the Ising spins. Si — ±1. In the end, the spherical model is formally the same as the Ising 
model with a global constraint instead of a local one: the sum of spins is constrained instead of each of them. At 
■ that time the saddle point method, needed in the solution, was not widely known, and here Berlin came in, leading 
"j--, , the celebrated jointpublication on the spherical model in 1952 Kac's personal reminiscence of this history is 
' presented in Ref. 0. 

^ , The spherical model for a ferromagnet has been considered in great detail. Actually, the paramagnetic to ferro- 
I magnetic transition is similar to an ideal Bose-Einstein condensation. Since the solution of the model is so simple 
■ O and explicit, the critical behavior can be solved exactly. Critical exponents and scaling functions can be derived. In 
^ ' particular, the model with short range interactions exhibits die = 2 as the lower critical dimension; for ci < 2 no stable 
Q ^ ferromagnetic phase occurs. Likewise, d^c = 4 is the upper critical dimension; for d > 4 critical exponents take their 
• • , mean-field values. These analytic results have been used to test approximations and general ideas of phase transitions 
^ ' for a wide range of interactions, short and long range. For a review on the classical spherical model see Ref. jsf. 

, As said, the spherical model was introduced for its mathematical simplicity. However, Stanley Q proved that the 
^ ■ free energy of a model of arbitrary spin dimension ly, incorporating thus the Ising model (for spin dimension ly ~ 1), 
' the X — y model (spin dimension v — 2) and Heisenberg model (i^ = 3), approaches that of the spherical model in the 
limit of infinite spin dimensionality v — > oo. Hence, it gives a geometrical interpretation to the spherical model. Since 
various critical properties where proven to be monotonic functions of the spin dimensionality v, the critical properties 
of the Heisenberg model appeared to be bounded on one side by those of the Ising model and on the other by those 
of the spherical model. 

The spherical model for antiferromagnets was thoroughly studied by Knops. The spherical constraint imposes 
< Sf >= 1 for ferromagnets. However, this does not work for antiferromagnets because of the lack of translational 
invariance. To recover this, he added a second constraint; more generally, one constraint has to be added for each 
translationally invariant set, which in the case of antiferromagnets means each of the two sublattices. Knops found that 
the two constraints reduce to a unique one provided the staggered external field vanishes. The fact that the spherical 
spins are scalars make it impossible to define an order parameter that can be identified with the spontaneous staggered 
magnetization. To solve that and get the proper order parameter Knops used a vector version of the spherical model 
0. He also generalized Stanley's arguments to non-translational interactions 
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The spherical model has also been applied to disordered systems. Though, in view of Knops' finding, perhaps an 
infinite number of spherical constraints should be used, typically no analog of the staggered external field is applied, 
and one may expect that all constraints collapse into a single one. Therefore spherical spin glass models may still 
give insight in the physics of the problem which would be more difhcult to study for e.g. Ising spins. In the case of 
pair couplings the exact solution exhibits no breaking of replica symmetry and the replica trick need not be used . 
The family of p-spin spin glasses (p-spin models) |S| has been shown to exhibit one step replica symmetry breaking 
by studying the spherical version. For spin glasses with random pair and quartet interactions {{p = 2} + {p = 4}, 
"p = 2 + 4"), one of us showed that an exact solution exists, exposing the full replica symmetry breaking scenario. 
The simplicity of spherical models thus may give insight in difficult problems for which otherwise no exact solution is 
available. For an early review on the use of the spherical model in disordered systems, see Ref. 0. 

So far the discussion has been classical. The classicality can be understood in particular because the entropy 
diverges at low temperature as InT, just as for a classical ideal gas. Different quantum versions of the spherical model 
have been proposed. Obermair studied surface effects in phase transitions [lOj. Identifying with a spin an operator 
Si, he postulated a momentum operator 11^ conjugate to it, [S'i,IIj] = ihSij. To get a spin dynamics, he added a 
kinetic term ^ ^ Il| to the Hamiltonian, but kept the constraint the same except for expressing it in operators 
as Tji{Sf) = N. In this case, the kinetic term may be understood as the kinetic energy of rigid rotors. The model 
remains exactly solvable. Many others have therefore used this quantization in the study of spin glasses systems 
with multispin random interactions (p-spin glasses) or the study of quantum phase transitions IT^. 

One of us presented in 1995 a different quantum approach to cure the low temperature behavior |l4j . In a Trotter 
approach to the partition sum, the first step is to take as spherical constraint: = Na/fi^ , where a is a constant 

that need not be unity, and Si is the complex parameter characterizing the coherent state associated with the bosonic 
annihilation operator = [Si/h+iYli]/\/2. Hence, in this approach the momentum appears in the constraint. Indeed, 
this constraint may also be written J^ii^f) + ^^(nf) = 2Na. As a second step, momenta dependent Hamiltonians 
were considered, by replacing JijSiSj JijS|i]j. Later the same formalism was applied to the p-spin glass 
model and was compared with its Ising counterpart. In spite of the simplicity of the system and its solubility, the 
resulting phase diagram shows very interesting critical phenomena. Since the momenta are present in the constraint 
one may also study situations where they does not appear in the Hamiltonian. 

Below we will consider two Hamiltonians with nearest neighbor ferromagnetic interactions. We will see that the 
presence or absence of momenta can change the symmetries of the action giving rise to different critical phenomena in 
the quantum region. The resulting action may be invariant under unitary transformations or under orthogonal ones, 
while in Obermair's approach only the latter is possible. In the last section we show that one of these spherical spin 
models relates to a quantum ferromagnet and the other to a quantum antiferromagnet. 

The two different quantum versions of the model have the same quantization rule [S'i,Hj] = ihSij. Both of them 
cure the problem of the entropy, it remains positive and, for temperature going to zero, goes to zero as a power law. 
Vojta, [13 following Stanley's arguments, found that Obermair's quantization gives a free energy that is identical 
to the large-n limit of the 0(n) nonlinear sigma model. Therefore it describes rotors instead of Heisenberg spins. 
Nieuwenhuizen, conversely, gave indications that his version had a behavior closer to Heisenberg spins, as having in 
the case of free spins a gap between the ground and the first excited state scaling with the field at small fields. 

The aim of the present paper is to point out the differences between these two models. We study two Hamiltonians 
using Nieuwenhuizen's spherical constraint. The first one was introduced in Ref. [T^ and we find it to be analogous to 
the large- A/" limit of a SU(A/') Heisenberg ferromagnet. In the second Hamiltonian studied, no momenta are present; 
momenta only appear in the formalism through the spherical constraint. In this case, the same critical phenomena as 
in the large- A/" limit of a SU(A/') Heisenberg antiferromagnet is found, which can be described by an 0(A/') nonlinear 
cr-model which, in turn, is analogous to Obermair's model. 

The paper is organized as follows: in section^ the classical spherical model is reviewed and the way to quantize it 
is discussed. The differences between the two quantum spherical constraints are pointed out. In section ITTll the path 
integral formalism is introduced to calculate the partition function of a quantum spherical model. In section Hvl this 
formalism is used to solve the thermodynamics of a ferromagnetic quantum spherical model with nearest neighbor 
interaction and the critical phenomena is studied in detail. At the finite temperature phase transition, the critical 
exponents remain the same as the classical ones. In section a Hamiltonian with the same couplings but without 
momenta is considered. The critical exponents are found to be the same as in Obermair's model. After that, in section 
I VII a generalization of the two types of Hamiltonians presented here is given and the limit of SU(A/') Heisenberg spins 
is argued to give the same critical behavior as this quantum version of spherical model. Finally some conclusions are 
drawn. 
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II. CLASSICAL SPHERICAL MODEL 

The spherical constraint was conceived as a relaxation of the Ising constraint. Indeed, Ising spins, Si = Si^z — if 
obviously satisfy it. Adjusting the coefficients from the original version it may be written as 



1 ^ 

\Y^S',=Na^ (1) 



i=l 



with a = ?i^/8 having dimension (Js)^. The Berlin-Kac spherical model is defined by the partition sum 



j DSe-P"5{\Y:^Sj-Na)^ J DS 



27ri 



where 



/oo 
dS, (3) 
-oo 



A. Vector spherical spins 

For vector spins the generalization of cq. |^ in the case of m spin dimensions reads 



2—1 a— 1 

It is worth mentioning that the spin dimensionality in eq. Q is not related to the approach of Stanley, who started 
with vector spins and ended up with scalar spherical spins. We only introduce vector spherical spins to avoid the 
restriction scalar spins have. We benefit from the fact that the vector character allows to study the behavior in a 
transverse field. A similar step allowed Knops to define a proper order parameter for the antiferromagnetic spherical 
model . 

B. Quantization 

It is natural to consider the Si analogous to position variables of harmonic oscillators. In quantum mechanics they 
become hermitian operators 5*^ with the dimension of h, 3s. The conjugate momentum operator Ilf is dimensionless 
and postulated to satisfy the commutation relation 

[St,fl';]^ihS,Ja,b (5) 
As for harmonic oscillators, this allows to define creation and annihilation operators 



^ = -^s^ " ^n- , s- = -^sf + -^fif (6) 



satisfying the commutation relation 



= (7) 

C. Spherical constraint on the length of the total spin 



There is some freedom to choose the spherical constraint, which amounts to describing different physical situations. 
The standard quantum constraint considered in literature is just the quantized version of the mean of eq. Q), 
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Constraint 1 



W<{Stf>^Nma (8) 



where < ... > denotes the quantum expectation vahie. Obermair took as the quantum Hamiltonian the classical H{S) 
with spins replaced by operators, and added the kinetic term that one expects for physical rotors, 

H{S,U) = ^gJ2fil + H{S) (9) 

i 

where is the rotor's moment of inertia. An effective Hamiltonian which includes the constraint can be derived 
with a Lagrange multiplier. One ends up with 

^tot = lgJ2llf + H{S) + ^,{t) [i ^(Sn' - Nma] (10) 



where /i is the Lagrange multiplier that enforces the constraint. In equilibrium its value is given by the equation of 
the spherical constraint ^ = 0. The dynamics is now fixed by the Heisenberg equations of motion, 

d 9° 

^ = t[H,,,{t),snt)]=gUf{t), (11) 
= ^[HUt),t^nt)] = -^-^<t)snt). (12) 

where the real parameter ^ has to be taken time-dependent in order to satisfy the "soft" constraint ((SJ at each instant. 
It is clear that a non-zero g is needed to get any spin dynamics. Combining the two equations one has 

g dt^ dSf 

Is worth remarking that no energy budget is involved in the spherical constraint, 

(iJtot) = {H). (14) 

D. Spherical constraint on the number of spin quanta 

In 1995 one of us has proposed a constraint that fixes the number of quanta |1J|]. In a path integral approach it 
was assumed that the c- numbers S", which characterize a coherent state, satisfy at each timestep 

Constraint 2' : ^ = Nm-^. (15) 

It is to be expected that this is equivalent to 

Constraint2: <^V^i < f^f >^ ^rn-^. (16) 

We shall show below that this is indeed the case, and the relation, a = a + h^, is derived in Eq. (|38|l . This constraint 
includes the momenta as can be seen by writing it in the form 



Constraint2: ^J2[{Sf^)+n'{uf))^Nmia+^) (17) 



I, a 



For a Hamiltonian i?(S, fl) that may, but need not, depend explicitly on the momenta, the effective spherical Hamil- 
tonian is 

H,,, = H{S,U) + \^,Y,[[S'^f + n\^i'tf] -Nm^,{a+^) (18) 
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Now, the situation where the Hamiltonian does not depend expHcitly on the momenta (no kinetic term), H{S, EE) 
H{S), still leads to sensible dynamics, since the constraint already depends on the momenta. Different constraints 
describe different physics. However, at high temperatures one expects the differences to become small. 

Eq. now brings Ilf — (dSf /dt)//i(i), eq. (|12|l remains the same. They may be combined together into a second 
order equation for the spin operators, 

In the remaining of this paper we will simplify the notation by taking units in which h = 1. 



E. Comparison of the two constraints 



The main difference between the two constraints is obviously the presence or absence of momenta. In the second 
case eq. the spherical constraint can carry all the dynamics of the model. On the contrary, using the first 

constraint eq. ||SJ|, a kinetic term, with an external parameter has to be added to the Hamiltonian .IQ] . This 
parameter determines the strength of quantum fluctuations; the classical model can be recovered for q = 0. This 
fact makes models with the first constraint describe quantum rotors, as was pointed out in Ref TS]. The first 
constraint, eq. ||SJ), brings actions which are invariant under orthogonal transformations. Conversely, using the second 
constraint, eq. the choice of Hamiltonian can bring symmetry under unitary transformations or orthogonal 

ones depending on the question whether the Hamiltonian contains momenta or not. Hamiltonians with unitary 
transformation symmetry yield free energies analogous to the large Af limit of the generalization of SU(2) Heisenberg 
spins to SU(A/'). Hamiltonians with orthogonal transformation symmetry share the critical phenomena with the large 
Af limit of 0{J\f) non- linear sigma model and describe therefore quantum rotors as occurs by using the first constraint, 
eq. ©. 

Each of the symmetries belong in different universality classes in the quantum regime, yet classical critical phenom- 
ena are always the same as in the classical model, consistent with the expectation that quantum effects do not lead to 
qualitative changes at finite temperatures. We will see that the dynamical critical exponent z is be different in both 
symmetries, causing the difference in critical exponents at the quantum critical point as was pointed out in Ref. [16| . 



III. PATH INTEGRALS 

In this section we explain, following Ref. [Tsl l , how to add the spherical constraint to a quantum Hamiltonian using 
the path integral formalism for models with the second constraint eq. H16|l . In second quantization the spins are given 
a bosonic algebra. In the path integral the boson coherent state representation is used for the spins (for a review of 
path integrals and coherent states, see e.g. 0| and for a complete study of coherent states see e.g. |l8|'l. 

A. Bosonic coherent state representation for a single oscillator 

Fock space is the Hilbert space of states labeled by the number of oscillator quanta. Coherent states are defined as 
the eigenstates of the anihilator operator d. Then it can be proved that for a system with many particles 

|^).e5:=^-i|0>=n(E^^^| 10) (20) 

a [ n„ J 

is a coherent state, where | ) is the vacuum representation in Fock's space, and a stands for each state for any particle 
of the system. Indeed, because of the identity daial^)^" = 'rea(al)"°~^ + (aJj)"°aQ, it holds that Ua \ (p) = (j)a \ (j))- 
The scalar product of two coherent states gives 



(</,|0') = eS=^^*^ (21) 
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A crucial property of the coherent states is that they form an overcomplete set of states. Any vector in Fock space 
can then be expanded in terms of coherent states. This is expressed by the closure relation [13 

fY^MMmM^.^^r.^^^^^^^^ =1 (22) 

''a 

where the measure in the integral comes from gaussian integration with complex variables (5 stands for imaginary 
part and 3? for real part) and the exponential term is due to the fact that coherent states are not normalized. Let us 
check that eq. (|22|l is indeed a representation of the identity of Fock space. We insert it in the left hand side of eq. 
(PT|l and we get 



( </. I 1 I 0' ) = 1 n '^^^^"^^^^^^"^ e- S = ^= ( I ) ( I ) = 

(23) 

a 

which indeed is the right hand side of eq. 121|) . 

The partition function of any quantum system Z = tr e^'^^'-" can be computed by the Trotter approach. The 
exponential has the same form as a time evolution operator in imaginary time. Thus it is possible to create a path 
integral over closed paths. The procedure is to split the exponential in a product of M equal terms. Between each 
pair of them a representation of the identity, eq. H22fl . is inserted. The partition sum then has the following shape 




where e = P/M and each 1 is an identity operator. Each of these identities is given an index; they represent the 
steps the system passes through in a discretized path. By using the identity defined in Eq. I|22|l the following matrix 
element is needed: 

{cl^,\e-^"^''''U^j-i) (25) 
Provided the Hamiltonian is normal ordered, the outcome is [T^ 

( 0, I e-^^('^'^'^) I ) « ( 0, I 1 - eH{a\a) \ <^,_i ) = e^*-^-^! " eff(^*,0,-i)) 

Correction terms can be neglected in the limit AI ^ 00 Each identity brings an integral at each time step. 

These integrals cover any path between its initial and its final state. The trace will finally tie the ends giving a closed 
path. The partition function finally reads 



[ i^(0*(T)^„(T))expj^dT 
-'0„(/3)=0c(O) [^^0 

where the subindex of the integral reflects the trace structure of the partition function since it gives a closed path 
integral; r stands for the imaginary time step, so 4>{t) — ipi', dr is the imaginary time difference between steps, so 
4>{t — dr) = 0i-i; and 



dm 

dr 



H{rir),^{T-dT)) 



(27) 



d(l){T) ^ (/)(t) - 0(r - dr) ^ (j>i - 
dr ^ dr " (3/M ^ ' 

Despite the fact that the nomenclature used in these formulas suggests a continuous time, it should always be 
understood as being discrete. The limit M ^ 00 should always be taken at the end of the calculations, otherwise 
some indeterminacies may arise. Continuous notation is used nevertheless because it is more compact. 
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B. Coherent state representation for spherical spins 

We can deal with spherical spins using almost the same approach. The operator a,; is identified with "Ef, where the 
index a denotes the spin vector direction, and the corresponding fields (pi are denoted as Ef. We remind that the 
spherical constraint we use is the one defined in eq. H16|) . 

In order to impose this constraint in the path integral formalism, the identity definition eq. H22() is modified to 
adopt to the spherical case, in a way inspired by Ref. 15]: one restricts the path integral to states which exactly 
satisfy the constraint by employing the truncated identity 



1 - l.p.e.w ^CJH dm'^)dmr) ^-i:'.i: | s )( S | <5(n - Nma) (29) 

ia 

where the number operator n, 



h = J2±f±l (30) 



counts the total number of spin quanta. We insert 



, 27r ./_™ 27ri 



— zoo 



where = z/2 is imaginary. (Strictly speaking, we should insert a Kronecker-(5 function, rather than the Dirac-(5, but 
for large N this amounts to the same.) Repeating the same procedure with this new identity we get 



Z= / D/iDS*DSexp(-A), (32) 

^S(/3) = S(0) 



with the action 



r=0 

and integration measures defined as 



'"s*(r) • + M(r)(S*(r) • S(t - dr) - Nma) + i?(S*(T), S(r - dr)) 

dr 



(33) 



/ Ds*£)s = n / / 



d5(E°(r))dSR(S°(r)) ^^^^^ 



ioQ 



/--n/::i^ 

so /i(r) is the Lagrange multiplier introduced to impose the spherical constraint and the prefactor C^*^-* is added to 
ensure, if needed, a proper normalization. Details on this factor were given in Ref. [Tsl l . 

It should be noted that the particle number operator in the definition of the identity eq. H29|l will be surrounded, 
as is the case for the Hamiltonian, by spin operators on different timesteps; therefore its creation and annihilation 
operators will also be projected on different timesteps, I]^*(r)I]f(T — dr). In Refs. 0,0] the spherical 

constraint was slightly different from the one presented here. The proposal was to take the constraint not in terms 
of the particle number operator but in terms of its generating variables (which are c- numbers), at every imaginary 
timestep, 



N 



^*-^ = Z.Z.^rWS«(r)=7Vma (35) 



so they acquire the same time-index. The two actions then differ only in the timestep projection of the spherical 
constraint, so with this constraint one obtained S*(t) • S(t) rather than S*(t) • S(r — dr). The difference that this 
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brings can be seen as follows. Starting from eq. H32|l we want to have two operators projected at the same time. To 
achieve this, it turns out that we must exchange the order of the operators, S^S = SS^^ — 1. The term EE^ can be 
projected at a single time as one can see following eq. (|26|l. For a single component spin at timestep j the relevant 
matrix element is 

(E,+i|e-^V^'^^'^e-^^|E,_i) - (S,+i|e-^^(l + e/i-e/iESt)e-i^|S^_i) 
= (l + eM)(E,+i|e-^^|E,>(E,|e-^|E,_i)-eM(S,+i|e^'^E|E,)(E,|Ete-^|E,_i^ 
= (l + 6M-eMS*E,)(E,+i|e-^^|E,)(E,|e-^^|E,-i)=e^^-^^^^'^^(E,+i|e-^^|^ (36) 

Thus a factor e'^~'^^> i"^) comes for each spin operator E^ at each timestep. This leads to the spherical constraint 
ls,Her^cal ^ / R ^^i^l^!^ e " ^ ' ' ^ | S ) ( S | c5(^ Er E," - Nm - Nma) (37) 

ia i,a 

that with definition cq. (|35|l should be compared to eq. I|29ll . In words, the spherical constraint can indeed be taken 
on the coherent state variables as in cq. 135|) or in Ref. 14], provided one makes the identification = cr + 1, or, 
restoring units, 

a = a + h^ (38) 

In Eq. (|48|1 we will verify that with this identification the two approaches indeed yield the same free energy. 

It is worth remarking that we imposed the spherical constraint strictly, no thermal average has been performed. 
In the following section /i will be integrated over by the method of steepest descends, a procedure that allows the 
particle number to fluctuate; therefore the satisfiability of the constraint remains only in average. 



IV. FERROMAGNETIC HAMILTONIANS WITH CREATION AND ANNIHILATION OPERATORS 

Using the formalism described in IIIII we can study the Hamiltonian 



V2 



i^jy(s,s,+n,n,)- Vr„;S,: 



(39) 



where in the second equality we inserted eq. ® . The iS'^IIj cancelled since we assumed symmetric couplings, Jy- = Jji. 
Obviously, the momentum operators do occur in this expression. The couplings Jy can in principle express any kind 
of interaction, ferromagnetic, antiferromagnetic, spin glass... The Ti represent an external field, that can be constant, 
variable, random... Later on, we will focus on ferromagnetic couplings in the presence of constant magnetic field. 
This Hamiltonian without the external magnetic field is symmetric under unitary transformations, a fact that will 
determine the critical behavior. 

The first step to get the partition function is to diagonalize the couplings, 



S,(r)=5]S,(r)e^ 

A 

S,(r)=5]S,(r)e,^ 



(40) 



where is the normalized eigenvector of the coupling matrix . 

Keeping in mind its ill-definedness, we may write the partition function sum as a continuum expression. 



/ 



DtDY,DJ:* exp 



^nr)—iy- + /i(r) (Er (r)E5;(r - dr) - Nma) 



- J.Er (r)E^(r - dr) -^T, (Er (r) + E^(r - dr)) 



(41) 
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In discrete notation, the action of eq. (|32(l reads 



a, A 



a, A 



-'Ai-AJ^Aj-l - 2^ i A ^ 



a, A 



a, A 



(42) 



where e = dr is the imaginary time step, j the time index and = is the field in the basis of eigenvectors 

of Jij . Collecting all terms we have 



A, a 



V 27ri 



^Aj-l^ 



V2 



where By = (Jy — (1 + e Ja — j+i! ^i^re the prime stands for the fact that (5^ m+i = ^ to the trace structure 

of the partition function. We can now integrate over the spins 



Z — \ Dfiexp 



^ mlndetBij + — ^^Bij"^ 



A, a 



(44) 



As usual, in thermodynamics, one-time quantities like /i(T) can be taken independent of r. We will employ this 
simplification throughout the rest of this paper. The determinant and the matrix inversion can then be performed 
[r^ . Integrating over /i by the saddle point method we obtain 



/3F = -ma/3M+^E|l^^(l 

A, a 



ax) 



2(1 -oa) 



(45) 



where a\ — 1 — e(/i — Ja)- Sending M — > oo we finally get 



2(/i - Ja) 



-Pfnn{a + -) + m / dJxp{Jx) ■{ In 



2sinh [ Jx) 



/3r- 



2(/i - Ja) 



(46) 



where in the last equality we have assumed that the couplings satisfy jfj^x'^x ^ ^- ^^'^ saddle point equation reads 



cr + 1 



iV ^ 1 1 - 6-/3 



r^A 



_e-'3(A'-JA) 2(^-Ja)2 



= / dJxp{Jx) 



1 _ g-/3(p-,/,) 2(/i-JA)2 



(47) 



The sums over the different eigenvalues of the coupling matrix have been changed into integrals. Each Ja has a 
weight in this integral given by p{ Jx)- The actual form for this weight function will depend on the type of couplings. 
A set of weight functions for ferromagnets in different cubic lattices can be found in Ref. ^ , and for spin glasses with 
long range interactions in Refs. 01 'H- 

In Ref. [Tsl l. where the spherical constraint used was the one in eq. (|35|) . the matrix B was different, namely 
By = (1 + ep,{je))6ij - (1 + eJA)'5(j+i- Then eq. ^ reads 



-Ppma 



N 



E{m 



2{p. - Ja) 



-Ppm{d - -) + m / dJxp{J\) <^ In 



2sinh ( ^(Ai - -^a) 



pr 



2{p - Jx) 



(48) 
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confirming that the already found shift a = a + 1, see eq. (|38|) . indeed brings the same value for the free energy. 
At large temperatures these equations reduce to 

(3F = -(3fima + m j dJxp{Jx) |ln/3(M ~ Jx) ~ ^{'^^-^Jx) } ^^^^ 

Apart from a factor two, these are exactly the equations of the classical spherical model, see e.g. Q. This factor 
two arises because the momenta double the degrees of freedom, see e.g. [l^. Near the phase transition they are 
already approximate, but the transition stays within the same universality class. 



A. Ferromagnetic couplings with transversal field in d dimensions 

In this section we will use the results given in the previous one for the concrete case of ferromagnetism couplings 
with uniform transversal field. The Hamiltonian in this case differs from the one before eq. (|39|l in the fact that 
the couplings only act in the z-direction while the external field only acts in the x-dircction (we restrict ourselves 
therefore to m = 2). The free energy reads 



[3F = -Pfi{2a 
and the saddle point equation 



W) 



■In 



2sinh(|(A.- J(k)) 



In 



2 sinh 



2 



2fi 



(51) 



where we have applied the changes Jx 



{2TrY 1 
J(k) and 



1 



g-/3(p-J(k)) I _ g-/3/. 



2/x2 



(52) 



dJxp{Jx) 



d<^k 



(53) 



We choose J(k) « Jq — J'|k|^ for |k| 0. In the case of short range couplings, for instance, one has x — 2 since 
J(k) — J2 Jcoski « J(0) — i J|kp. A long range couphng that decays as J{r) ~ 1/r^" at large r gives x — a — d. 

As in the theory of Bose-Einstein condensation, the saddle point equation fixes the dependence of ji on temperature. 
There should be a solution at any T. In order to have a real free energy, p cannot be smaller than the maximum value 
for J(k). Therefore, we should investigate the convergence of the integral in the limit fj, ^ Jq. If the integral diverges, 
(3 must go to infinity before fi reaches Jo in order to satisfy the saddle point equation, so there exists a fj, for all 
temperatures and no phase transition occurs. If the integral converges, however, there will be a range of temperatures 
in which the saddle point as it stands cannot hold. This indicates that we have overlooked a macroscopic occupation 
of the ground state, as occurs in Bose-Einstein condensation. The relevant integral behaves as 



d'^k 



1 



1 



dk¥^' 



(54) 



where VLd is the hypersurface of a sphere in d dimensions. At k = 0, this integral converges for d > x, hence there 
will be a phase transition for dimensions larger than x. 

At low temperatures, /i may get stuck at Jq and the saddle point equation as it is in eq. (|52|l is no longer valid. This 
is because, as in Bose-Einstein condensation calculations, the ground state is not properly included in the integral. It 
should be taken out of the sum before this one is converted to an integral. This causes a change in the free energy by 
a factor {jj, — Jo)g and the saddle point equation becomes 
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where q = jf {^k=o'^k=o) ground state occupation, q can be evaluated from the saddle point equation (/z — 

Jo)^/q — 0. Thus when n = Jq the occupation of the ground state can take non-zero values that can be determined 
using eq. (|^ . Hence the ground state occupation is macroscopic in the ordered phase. 

A transversal field will lower the transition temperature. Above a certain value Fc, the transition does not exist 
anymore, thus T = 0, F = Fc is a quantum critical point (for a complete study over quantum phase transitions, see 
e.g. 0]). We will now first study the classical critical point, where F = 0. 



1, Finite temperature phase transition 
For the dimensions where the phase transition exists, the critical temperature is found by solving the equation 



2(cr + l) = 



1 



1 



1 



-Pa Jo 



(56) 



The dependence of the chemical potential on the temperature near the transition is the first thing needed. To get 
it, we expand the saddle point equation around the critical point T = Tc + t, ^ = Jq + Sfi. The integral gives, up to 
first order in Sfi and r 



^ (27r)'i 1 - e-'3(A'-'^(k)) 



" d'^k 



1 



Jo - J(k) 



1 _ p-/3e(./0-J(k)) 



4T^ sinh' 



(^) 



— Sfi- 



1 



ATr sinh 



2 ( ./n-./(k) 

2Te 



(57) 



The coefficient of 5fi is an integral that diverges for d < 2x. This means that for these dimensions the leading term 
in the (5/i expansion of eq. (|57() has a power smaller than one. For dimensions d > 2x we will have (x t which will 
lead to the mean field exponents, 2x is therefore the upper critical dimension. To study the system near the critical 
point we subtract eq. I|56|l from the saddle point equation, a procedure that will cancel the zeroth order term in the 
expansion in t and Sfi, giving finally 



where 



T ~ ad=2xS^J.ln5fi, ad=2x = a 



ArtdT^n 



T w Sd>2xSfi, ad>2x = aTc 



(2^)'^J'^a;sin(i^^) 

{2ttYJ'^x 

d'^k 1 



(2-)%inh^ (:^) sinh^(^) 



for X < d < 2x 



for d = 2x 



for d> 2a 



d'^k Jo - J(k) 



Jo 



V--(2'^)'sinh2(i!i^) sinh2(^ 



is a finite, positive number. 

The internal energy of the system reads 



(58) 



(59) 



U = -/i(2(T + 1) 



d'^k At-J(k) 
(2^ 2 



coth 



- J(k)) 



— coth 



2^ 



(60) 



The specific heat close to the transition from the paramagnetic side can be written as 



for X < d < 2x 
for d > 2x 



(61) 
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where 



1 r A [Ai-J(k)] 



2 



-2fT- 1 



(2tt¥ 



- coth 



Ai- J(k) 



2T 



/i - J(k) 



4Tsinh^ 



2 I^ M-J(k) 
2T 



2 \2TJ 



4rsinh2 (^) 



r2 



where are the prefactors in eq. H58|) for the corresponding dimension. In the ordered phase fj, is stuck in its 
minimum value (fi = Jq) for any temperature. Hence, C = Co{fi = Jq) in the ordered phase. The critical exponent 
a is the expected one: a = ^jE^ for x < d < 2x, and the mean field value a = holds for d > 2x, which describes a 
jump in the specific heat. 

Adding a small longitudinal field h, the free energy reads 



PF = -/?^(2cr + 1) 



{2TrY 

and the saddle point equation becomes 



In 



2sinh( J(k)) 



In 



2 sinh 



2 



2/Lt 



2(m - ./o) 



(64) 



2(a + l) = 



{271^ 1 - e-/5(A'-./(k)) ^ 1 - 



2/i2 



(65) 



By differentiating the free energy with respect to h can be seen that the magnetization is = jj-rj^ ■ In the limit 



3 with 



/i ^ 0, it is proportional to the square root of the occupation of the ground state, since by comparing eq. 

eq. l(55|l one finds q = -^(E^J^qS^^q) = The factor i appears because is actually the real part of the spin field 
the one macroscopically occupied and a half term appears in the change eq. ©. From eq. H65|l . we can approach 
the transition by sending the longitudinal field to zero at the critical temperature. The saddle point equation now 
accounts for the dependence of the chemical potential on the field. The calculation is similar, yielding finally 



2rtjTcTT 



(2^)''J'''/"a;sin(^^) 
2VLdTc 




In 5^ 



(2^)' 2r.sinh2 



5^ 



3/2 



for X < d < 2a; 



for d — 2x 



for d> 2x 



(66) 



(^) 2r.sinh^(^) 

Therefore the critical exponent 5 is given hy 5 ~ for dimensions x < d < 2x and the mean field value (5 = 3 is 
recovered for d > 2x. From the magnetization, the susceptibility follows as x « Therefore we find 7 — for 
X < d < 2x and 7 = 1 for d > 2x. In the ordered phase, the expansion of the saddle point equation, eq. H65() . for T 
near the transition yields 



M 



2T^ 



^ d'^k Jo-J(k) 



Jo 



. (2-)%inh^ (irfs)) sinh^(^) 



(67) 



Therefore for all dimensions where the phase transition exists, one has (3 = 



13 



For other critical exponents the correlation function is needed. It can be computed adding the right source term to 
the Hamiltonian, ^xi9>^('''q)^\. ~^ 9\i'''r)'^\) a-nd differentiating 



92 



d9l{Tq)dgy{Tr) Za 



(68) 



5* =5=0 



where T stands for the time ordered product. Z{g*,g) is the partition function of the Hamiltonian including the 
source terms and Zq is the partition function without them. This procedure is carefully explained in Ref. (l7l | giving 
the result 



G(fc, t) = G(fc, r|fc, 0) = e^^^-'^C^" {9{t ~ r]){l + Uk) + e{-T + r;)nfc} (69) 

where is a Heaviside step function and 77 is a positive infinitesimal that indicates that the second term is the relevant 
one at T = 0. Furthermore, 



nk = -5 (70) 

is the boson occupation probability, with uj — (/i— J(k)). Fourier transforming this last result to Matsubara frequencies 
we get 

So when we approach the critical point, we can see from this equation that oc S/i, then using eq. H58(l we find 
that = for dimensions x < d < 2x and v = ^ for d > 2x. rj — 2 — x due to the fact that the couplings depend 
on , and z — x because in the denominator lo appears as a linear term. Both rj and z are valid for any dimension. 

This finally gives all the critical exponents of this finite temperature phase transition, which are exactly the same as 
in the classical model. This is expected from renormalization group arguments [l6|]. The critical behavior is controlled 
by a classical fixed point, therefore quantum dynamics does not play a qualitatively new role. Hence, the results are 
the same as in the classical spherical model Q or other models with different quantum dynamics considered at finite 
temperatures | il3j . 



2. T=0 quantum phase transition 

In this section we analyze the behavior of the system at T = 0. As it can be seen from eq. (|52|) . when the transversal 
field increases, the temperature of the transition decreases till it reaches zero. This defines a quantum critical point 
Tc = at r = Fc. In order to study it, an analogous procedure as before should be followed. At T = everything 
happens to be rather simple. The free energy reduces to 

F = -2a^i - ^ (72) 

The saddle point equation turns out to be 



in the paramagnetic phase, 

in the ferromagnetic phase. (73) 

where q = ^M^ is the occupation of the ground state for small transversal fields. Since the temperature vanishes, 
quantum fluctuations, controlled by F, give rise to the phase transition. Therefore, the parameter that should be 
used to control the transition is the transversal field and not the temperature. Then the proper analog of the specific 



2(7 = 



2^*2 



2(7 



2J2 
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heat will be proportional to the second derivative of the free energy with respect to the source of fluctuations, the 
transversal field. 

^^r = ^--- + ^^ (74) 

As before, we must know the dependence of Sfj, (fi = Jq + 5^) on the distance to the critical point {5T) in the 
paramagnetic phase. From eqs. ()74I73|) . it can be seen that the analog of the specific heat has a jump discontinuity, 
implying a — Q. 

Cr — in the paramagnetic phase, 
Cr = —f- in the ferromagnetic phase. (75) 

where the minus sign comes from the fact that the T = free energy, eq. (|72() . is negative. Adding a small longitudinal 
field as before, we find the critical exponent (5 = 3, since oc ^ and J/i (x at Tc- For the susceptibility, we find 
7 = 1 since x oc 6ijl~^ oc (5r~^. Subtracting the saddle point equation near the transition in the ferromagnetic phase 
from the one at the transition, we get q = ^ „ , which in the lowest order gives Ml cx 5T and therefore (3 = h. Eq. 

(|71|l can be used here once it is transformed to real frequencies, iw„ — uj + irj. Then since S/i oc ST we find v = 
rj = 2 — X and z = x. Hence we find always the mean field values. This occurs because the quantum critical point of 
a d-dimensional model shares the critical exponents of a classical critical point of a (d+z)-dimensional model, as it 
was shown by general renormalization group arguments |l(i| |. Then, since the transition exists for d > x, the quantum 
critical point behaves as the classical one for d > x + z = 2x. 2x is just the classical upper critical dimension in this 
model, above which one has the mean field values. 



V. HAMILTONIANS INVOLVING SPINS BUT NOT THEIR MOMENTA 

In this section we are going to extend the analysis of section Hvl to a Hamiltonian which only depends on the spin 
operators S and not on the momenta 11. When going from the classical to the quantum model, we have to keep in 
mind that the Hamiltonian must be Hermitian. To be precise, the Hamiltonian we will deal with is 

iJ = -i^J,,S,S,-^r,S. (76) 

ij i 

with real valued Jij and and where S ~ (S^ + 1]) / \/2 is the real part of the former spin field. Hence, the Hamiltonian 
does not involve momenta, but the spherical constraint does, see eq. ()16|l . This changes the symmetry of the problem 
from invariance under unitary transformations to orthogonal ones. In terms of boson creation and annihilation 
operators the coupling term for symmetric interactions, Jij = Jji, is proportional to Jij(2S|Sj + S^Sj + S|sj), 
where we can notice the symmetry of the problem. We will see that this model reproduces the 0{Af) quantum rotor 
model. 

We can get the partition function in many ways. A similar procedure using discrete imaginary time path integrals 
can be done as before. This gives us many problems due to the fact that creation and annihilation operators are 
projected on different time steps which is a lengthy and tedious procedure. However, the form of the Hamiltonian 
makes it suitable to apply a Bogoliubov transformation (for details see e.g. H^)- Due to that, we get the same 
Hamiltonian as before but with different coefficients. In order to do so, we must add the spherical constraint directly 
to the Hamiltonian. The procedure is as follows: first the couplings matrix is diagonalized by inverting the lattice as 
done before in eq. H4U|) and then the 5"s are shifted to absorb the field term. This finally gives 



Performing the Bogoliubov transformation it turns into 



Nmiia, (78) 
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which is a Hamiltonian analogous to eq. H39|) . So it can be diagonaHzed as explained, giving finally 



(3F = -(3^im{(7 + -) + m I dJxp{J\) \ In 



2 sinh ( - y^Ji(Jj~7x) 



/3r 



2(/i - Ja) 



(79) 



where we have put back the factor m standing for the number of components of the vector spin. The saddle point 
equation is obtained as 



a + 1=1 dJxpiJ,) I J^^^ coth f VMa^-J.) + 



4^ pip - Ja) 
At large temperatures these equations reduce to 



2ip-Jxy 



(80) 



PF = -ppm,{a + dJxp{Jx) (^\n Py^ p{p - Jx) - ^^^J^ 



(81) 



''^l = I''''^''^{fp^2U;hx}' 2ip-Jx? 



(82) 



These equations are very similar to the standard ones of the classical spherical model (up to a factor 2), see eq. 
(I49|l . but they are only identical where they should be, namely at large T, where also p ^ T \s very large, see also 
Ref. ^. 

A. Ferromagnetic couplings in the presence of a transversal field 

Analyzing the phase transition of the Hamiltonian that does not depend on the momenta is analogous to the 
previous case. We begin again by choosing the coupling term in the z direction and the external field in the x and 
we assume ferromagnetic couplings. The saddle point equation gives a phase transition via a macroscopic occupation 
of the ground state, which in the present case is a bit more complicated. The critical exponents are different, due to 
the fact that the symmetries of the system have changed. The free energy reads 



(3F = -(3p{2(7 + 1) 
and the saddle point equation 



{2nY 



In 



2smh{'-^pip~J(k)) 



In 



2sinh(^ 



2p 



- ^ (83) 



4cr + 2 = 



d'^k 2p - J(k) 



i^nr 2^p{p-J{k)) 
We now analyze this model in detail. 



coth 



Iv'M/^- j(k)) 



coth 



(3p 



//2 



(84) 



1. Finite temperature phase transition 

Following the same procedure as before we can find that the transition exists for d > x and that the upper critical 
dimension is d = 2x. The critical temperature is the solution of 



4(7 + 2 = 



d'^k 2 Jo - J(k) 



coth 



^V^o(Jo- J(k)) 



coth 



/3Jo 



(85) 



The dependence of the chemical potential in the temperature near the classical critical point reads 
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where 



ad>2x = OL 



ad<2x = a 



ad=2x = a 



2ndTj^n 



2ndT^ 

d 7/2 „ 



for X < d < 2x 

for d — 2x 
for d > 2x 



d'^k 2n - J(k) 



d 



dfi [J (27r)'i 2 - J(k)) 



; coth 



y/flifl - J(k)) 



coth 



d-'k 2 Jo - J(k) 



Jo 



- (2'^)' 2 sinh^ TMgZM^ sinh^ (^) 
The internal energy of the system reads 



(86) 



(87) 



rr (0 r "^'^ VMM- J(k)) ^, 



^VMa*- J(k)) 



^cothf^V?^ 



^2 
271 



The specific heat has the same expression as in eq. (|61|) where now correspond to the prefactors of eq. (|V A l\i 
and with coefficients 



Ci 



Co 



1 r d'^k /^(^-j(k)) 



4T2 {2nY ^.^^2 1^ V^(^M^ sinh^ (^) 



VMM- j(k)) ^2^'^^'' 



2n - J(k) 



STsinh^ 



2 / VM(M-J(k)) 



2T 



+ - coth f — ) - 
2 \2TJ 



4rsinh2 



(89) 



(90) 



This is anafogous to the previous model and gives the same exponent, a = for x < d < 2x and a = for 

d > 2x. Adding a small magnetic field longitudinal to the couplings, the free energy becomes 



PF = -/3M(2a + 1) + I ^ 



In 



2sinh( |Vm(m- J(k)) 



In 



2 sinh 



(3^i 



2n 2{fi - Jo) 



(91) 



therefore the magnetization is AIz — ^^j^ , which is as before the square root of the occupation of the ground state 
q — jj- < S'z(|k| = 0)^ >= AJ|. The saddle point equation is now 



2(2cr + 1) = 



TT ,d 



d'^k 2)1 - J(k) 



coth |^Vm(a*- J(k))| + coth (^^^ + + 



(m - Joy 



(92) 



With all these and following the algebra of the previous section one finds the same critical exponents for the 
magnetization for the same dimensions since we are in the classical critical point. 
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The time ordered correlation function (TS'^(t)S'Ij.(O)) differs from the previous one, eq. H69|l since in this case 
the S are not the variables that diagonalize the Hamiltonian in eq. (|78(l . We must write S in terms of a and then 
compute the correlations. This brings 



Gik,T\k,0) 



J(k) 



where 



VAi(/i- J(k)) 



Uk cosh 



which in frequency space reads 



r|^M(M-J(k)) 



(93) 



(94) 



G{k,iujn) 



-J(k) 



(95) 



2K-^(^- J(k))] 

When approaching the critical point we find that oc Sfi as before and we get the same value v = for 
dimensions x < d < 2x and = ^ for d > 2x. Since couplings appear in the same way as before we also get the same 
value, T] = 2 — X for all dimensions.. The difference appears in the dynamical critical exponent. Here cun appears 
squared, therefore z = x/2. Here we see how the model reproduces the critical exponents of the rotor model as in 
Ref. |1^] bringing thus a different behavior at the T = quantum critical point from the model of section llV A 21 



2. T=0 quantum phase transition 



In this case the T = phase transition is more interesting due to the fact that the dynamical critical exponent 
z = ^ is smaller than z = a; of the previous section. The free energy reads 



F= -^(2a+l) 



2 2 2^ 



and the saddle point is set by 



4(7 + 1 = 



{2t:Y 2 
d'^k 2^ - J(k) 



{^^Y 2v/^4;r~7(k)y 



(96) 



(97) 



We find that the transition exists for dimensions larger than c? > f and c? = ^ is the upper critical dimension. The 
chemical potential depends on the source of fluctuations iJF = F — Fc as 



5T K ad^3j./2(5/i In (5^, 
8V K. ad>3x/2SfJ-, 

0'd=3x/2 



^d<3x/2 



2{2ttYJ'- {2d-x)xy/n 



X , 3x 
ior - <d < — 
2 2 



i'2nr SxJ'i 



for d ~ — 
2 

tor a > — 
2 



(98) 



"^0 



2Vi 



d'^k 



2 Jo - J(k) 



•^0 7 {2'kY \^Jo(Jo- J(k)) 4(Jo(Jo- J(k)))^ j. 

where the F's on the right hand side of the first equality are Euler's Gamma functions. The specific heat (see eq. H74|l 
) coming from the disordered region will behave as 
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where is the prefactor in eq. (|V A2|) for the proper dimension. Coming from the ordered region, conversely, 
Cr « Therefore a = for § < d < 4f and a = for d > 

jQ za — X A J. z 

The dependence of a small longitudinal field on the chemical potential, in case the transversal field is at its critical 
value, reads 



h : 
h: 



ad=sjE. In (5^ 



for-<d<- 

for d = — 
2 



2r| 



for d > 



3x 



(100) 



From these equations and the ones for the magnetization and the susceptibility we can find that S ~ '^n^^iH and 



2d~x 



1 



for I < 0? < ^ , while (5 = 3 and 7 = 1 for o? > ^ . As before /3 = ^ for every dimension. For the correlation 
function the calculation is the same as in the finite temperature case, projected into real time, the exponents are 



2d-a 



for dimensions | < d < ^ and u = - above the critical dimension, r] = 2 — x and 



for all dimensions. 



VI. GENERALIZATION AND MAPPING FROM HEISENBERG SPINS 

In this section we generalize the two preceding Hamiltonians and we map the Heisenberg model onto the spherical 
model. In a more compact way, we can write the former Hamiltonians in absence of external field as 



B;, 



(101) 



If the matrices Aij and Bij can be diagonalized simultaneously, the techniques from previous sections can be used. 
The free energy reads 



/3F = -/3^m(a + i) + ^^ 

A 



PA, 



ln^2smh['-^/{pi-A^)^-Bl 



(102) 



where fi satisfies the saddle point equation 



a + - = - y ^'^^ coth/^./(^-A,)2-i? 

2 N ^ 2^{^i - Aa)2 - Bi 12V^^ ^' 



I 2 



(103) 



The coefficient Bij in eq. (|101|l is the responsible for a change in the symmetries of the problem. If Bij is zero, the 
action is symmetric under unitary transformations while if it is non-zero the symmetry is reduced to orthogonal. 

The mapping from Heisenberg spins comes as follows. The Hamiltonian can be written in terms of Schwinger bosons 
[20|. The Schwinger boson transformation for SU(2) spins reads 



= a\a2, S ~ aiCj, S^' — —{a\ai — 0202) 



(104) 



This can be generalized to SU(A/') spins and expand around the large- A/" limit |2J. In a path integral formalism 
for ferromagnetic interactions 



H 



2 E Jij^i ■ 



'2^ E J^^^. 



(105) 



where i,j represent lattice sites and m, n represent the boson flavor. The Hilbert space spanned by Schwinger bosons 
is much larger than the one given by Heisenberg spins. The constraint needed to restrict it to the physical Hilbert 
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space is that the number of Schwinger bosons at each site has to be kept fixed J2m '^m = J^S. This is inserted into 
the formahsm in the same way as we have done it with the spherical constraint, a Lagrange multipher Hilr) appears. 
The biquadratic terms can be decoupled by a Hubbard-Stratonovich transformation (see e.g. In the case of a 

ferromagnet, the transformation at each time step and for each flavor in the path integral reads 



exp ^ Jzja*aia*aj \ J Y[ dQ^j exp < ^ ^ QtjJtjQji ~ ^Yl Q^jJii'^j^^ } (106) 



where a field Qijir) has been generated. In the mean field approximation, one puts Qijir) — Q and fJ-ii^r) — fi. 
Hence, one gets up to a non interesting constant 

i , m ij, m ij 

where we have aheady added the Schwinger boson constraint. The free energy per particle reads 

(3F = ^Y1 1^(1 - e-^(^-'3-^('')) + J(k - 0) - fiSfjM (108) 

k 

and the saddle point equations are 



^E"k = 5 (109) 

k 

^ E -^(i^K = Q J(k = 0) (110) 

k 

where Uk is the boson occupation number eq. H70|l with ui — ^ — (3J(k). Subtracting the two saddle point equations 
we can see that for large S and small T we can approximate Q ^ S recovering then the spherical model eq. H46I47|I 
for zero external field or eq. (|102I103() for Bij = 0. From this approach we thus see that the free energy of a SU(A/') 
Heisenberg ferromagnet for large Af is formally the same as the quantum spherical model proposed in eq. 1)3911 in the 
thermodynamic limit, so when the radius of the hypersphere that defines the model {N in eq. H16|l ) is also very large. 
Thus the large Af limit is somehow analogous to Staley's large spin dimensionality limit. 

In the case of an SU(A/') antiferromagnet the procedure is more or less the same but the symmetries are different. 
The lattice is divided in two sublattices A, B. In one of the sublattices a spin rotation is performed that allows us to 
write the Hamiltonian in the form [2ll | 

= 2 E "^ii^i ' ^3 ~* ~ 2iV ^ '^ii^*im^*im'^jn'^jn (HI) 
ij ij.mn 

Performing a Hubbard-Stratonovich transformation as before, the Hamiltonian with the Schwinger boson constraint 
in the mean field approximation finally reads 



HiT'^'W = E - I E J^Mm4m + «™a.>n) + ^ E ^^^^ " A^^^/X (112) 



It is important to stress that here the SU(A/') symmetry has been reduced to a residual 0(Af). The free energy per 
particle reads 



/3F=:^El^(2sinh 



TV 

k 



/3AA(5 + lV + ^:^J(k = 0) (113) 



20 



and the saddle point equations read read 

1 J^(k)Q 

where n/j is eq. (|70ll for uj = \/ ~ J^(k). Subtracting the first equation times /i from the second times Q we get 

^ E v//^' - O'-^'Ck) ("k - = M + ^) - J(k = 0) (116) 

The first term is proportional to T, so for very small temperatures and very large S', near the transition where 
^ K, QJCk — 0), we can approximate Q ~ 5 + i. Then eqs. (|113I114|) are analogous to eqs. (|102I103|) for Aij = 0. 
This will have the same critical behavior as the model in section due to the fact that it comes from the term 
coth[\//L( — J(k)]/ ^/Ji^^~J(k) which also appears here due to the equality 2nk; + 1 = coth[A/ (/x + QJ(k))(^ — QJ{k.))]. 



nk+2 



QJ(k = 0) 



(114) 
(115) 



VII. CONCLUSION 



In this paper we have explained a way of working with quantum spherical spin models using path integrals and 
coherent states. Some examples of the use of this formalism are given, eqs. (|39I76|I . and their critical phenomena 
are studied. We propose a comparison with SU(A/') Heisenberg models that gives a geometrical interpretation to the 
quantum spherical spins. The spherical constraint we use, fixes the number of spin quanta E, eq. (|l()|l : in other 
words, it fixes both the average length square of the spin operator, S*^, and the one of its conjugate momentum, 11^. 
The usual version of the quantum spherical model, on the contrary, involves only the spin part S. The presence of 
momenta in the spherical constraint allows the Hamiltonian to have no kinetic term, since it can be induced by the 
constraint, a fact that can change the symmetries of the problem, and due to that, the critical behavior. 

The Hamiltonian in eq. (|39l) yields an action invariant under unitary transformations. It brings formally the 
same free energy as a SU(A/') Heisenberg ferromagnet in the limit of large Af. The other Hamiltonian studied, eq. 
()76|l . brings an action invariant under orthogonal transformations; it gives the same critical behavior as an SU(A/') 
Heisenberg antiferromagnet in the limit of large A^, which is, in its turn, analogous to an 0(Af) nonlinear cr- model 
or quantum rotor model |l3l [T9l | . The main difference between these models lies in the dynamical critical exponent 
z which brings a different behavior at the quantum critical point. Classical critical phenomena are, as expected, the 
same in both models and equal to those of the classical spherical model. 

In the formulation of the model, the strict spherical constraint has been used where fluctuations on the particle 
number are not allowed. The constraint is added to the action via a Lagrange multiplier. The strict approach has 
to be abandoned when we integrate this Lagrange multiplier using the saddle point approximation. In this step, 
we automatically allow fluctuations on the particle number and therefore the constraint ends being satisfied only in 
average. These effects are immaterial in the considered thermodynamic limit, but do enter finite size corrections. 

The analogy of the two Hamiltonians studied here with Heisenberg models in the large spin dimensionality limit 
has a drawback. Both models have different coupling to the external field. In spherical models it comes in linearly, 
as a source term. No analog has been found for this in the large spin dimensionality limit of the Heisenberg model 
where each spin contribution brings a bilinear term in Schwinger bosons. 

Another approach could have been to start directly from the SU(A/') Heisenberg model and to do the already stated 
large Af limit to get to a solvable model. In order to have a transversal field that competes with the ordering of the 
interacting spins one could introduce anisotropy in the model. A study of this type has been done for 2 dimensions 
by Timm et al. in terms of Schwinger bosons and in terms of Holstein-Primakoff bosons by Kaganov et al. |2^ 
for any dimension. The anisotropy term brings a residual spin symmetry describing Ising or XY spins. The phase 
transition depends on the type of this residual symmetry; an additional transversal field decreases the transition 
temperature towards zero giving a quantum critical point, result qualitatively reproduced by our model. 

In spite of the lack of direct interpretation of the source term in the mapping from Heisenberg spins, the phase 
diagram follows the expected behavior for a spin model with an external transversal field. The critical exponents 
for the classical and the quantum model are the ones expected by renormalization group arguments. The quantum 
critical point behaves as the classical one for dimensions Dquant — ddass + z where z is the dynamical critical 
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exponent. For instance, in the model with unitary symmetry eq. H39|) . z is just the difference between the lower and 
the upper critical dimension, thus the quantum critical point has the same mean field exponents as the classical finite 
temperature critical point. 

Many other models have a clear analog with this one. Sachdcv and Bhatt j2J| represented pairs of spins in a square 
lattice with a bond representation; they form either a singlet or a triplet. These elements can be written down in 
terms of the canonical "Schwinger boson" representation of the generators of SU (2) ® SU{2) — 50(4). Since a couple 
of spins either form a singlet or a triplet, a constraint must be added s^s + ^^ ^I^q = 1, where s represents the singlet 
annihilation operator, and ta represents a triplet annihilation operator in the a direction. Sachdev and Bhatt study 
using this formalism systems with interactions up to third nearest neighbors. They make the further assumption that 
the singlet part condenses and replace the s operator by its mean field value (s) — s, and solve the rest for the triplets. 
The final Hamiltonian is very close to our eq. H76|) . or, better, the generalization of our model eq. (|HJ1|I with the 
proper couplings. A minor difference is the role played by the non-constant mean value of the singlet part. 

An interesting line for future research would be to expand this version of quantum spherical spin models to different 
types of interactions and fields. Randomness is easily added in the model. The dynamics could also be studied for 
the constraint described here. Another line would be to study in this approximation the Heisenberg SU(A/') model for 
any dimension with anisotropy in one dimension and transversal field, giving special attention to critical phenomena 
and compare it to the Holstein-Primakoff approximation of Ref. . 
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